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Abstract We have studied ‘*He confined in a 95% porosity silica aerogel in the 
vicinity of the bulk liquid gas critical point. Both thermodynamic measurements 
and light scattering experiments were performed to probe the effect of a quenched 
disorder on the liquid gas transition, in relation with the Random Field Ising Model 
(RFIM). We find that the hysteresis between condensation and evaporation present 
at lower temperatures disappears at a temperature 7),/, between 25 and 30 mK be¬ 
low the critical point. Slow relaxations are observed for temperatures slightly be¬ 
low Tch, indicating that some energy barriers, but not all, can be overcome. Above 
Tch, no density step is observed along the (reversible) isotherms, showing that 
the critical behavior of the equilibrium phase transition in presence of disorder, 
if it exists, is shifted to smaller temperatures, where it cannot be observed due to 
the impossibility to reach equilibrium. Above Tch, light scattering exhibits a weak 
maximum close to the pressure where the isotherm slope is maximal. This behav¬ 
ior can be accounted for by a simple model incorporating the compression of ^He 
close to the silica strands. 

PACS numbers: 67.25.dr, 05.70.Jk, 64.60.My, 68.03.Fg 


1 Introduction 

Despite intense theoretical and experimental efforts, the effect of quenched dis¬ 
order on phase transitions remains not fully understood. While this question has 
originally arisen in the context of magnetism, it has been extended to fluids con¬ 
fined in disordered porous materials, in particular tenuous silica gels or aero¬ 
gels which offer an unique experimental realization of a self-sustained, dilute, 
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quenched disorder. Focusing on the Ising universality class, Brochard and de 
Gennes have suggested that the demixtion of a binary liquid conhned in a gel 
should be an experimental realization of the Random Field Ising Model (RFIM) 
and of its critical behavior ^ at least close enough to the bulk critical temperature 
Tc for the fluid correlation length to be larger than the gel correlation length. In 
the same spirit, a number of studies have concerned the liquid-gas transition of 
^Fle confined in aerogels. Early experiments by Wong and Chan— have been inter¬ 
preted in terms of a genuine first-order equilibrium phase transition existing in a 
wide temperature range, with a shifted critical point and a modified critical expo¬ 
nent for the order parameter (the density difference between liquid and gas), due to 
the disorder. However, latter experiments on '^He in aerogel s have revealed a 
strong hysteresis between condensation and evaporation, persisting up to close to 
the critical temperature* , indicating that the system cannot reach equilibrium. A 
remarkable feature is that, at low temperatures or large porosities, the adsorption 
isotherm presents a quasi-plateau of pressure, just like for a regular first-order, 
equilibrium, transitio n In fact, this behavior turns out to be consistent— with 
the prediction originally made for the zero temperature RFIM—, and generalized 
to the case of the liquid-gas transition at finite temperatur e^ — '^ ^ , of a disorder- 
driven, non equilibrium, phase transition, where condensation occurs by athermal 
avalanches, whose size diverges below a critical temperature (or disorder). 

These results leave open the possibility that, very close to the energy bar¬ 
riers for condensation become small enough for the system to be described by the 
equilibrium RFIM, characterized by a true phase transition below a critical tem¬ 
perature T*, and non bulk critical exponents at T*—. Experiments on N 2 — and 
isobutyric acid-water'^ confined in (aero)gels have been interpreted as showing 
such an equilibrium behavior, while the observation of slow dynamics above T* 
in these experiments fitted with the initial theoretical suggestion that, above T*, 
the Ising REIM model could present an equilibrium glassy phased. However, it 
has been recently rigorously shown that only an out of equilibrium glassy phase 
is possiblei^ above T*, and it is not clear how the experiments fit this picture. 

In contrast, two other experiments And no evidence for any equilibrium critical 
behavior. Eor ^He confined in aerogels, isotherms measurements very close to 7)., 
in the region where the hysteresis disappears, show no true pressure plateau, i.e. no 
indication of an equilibrium phase transition—. On the other hand, neutron scatter¬ 
ing measurements performed on CO 2 confined in a silica aerogel close to the bulk 
critical poin t show that the fluid correlation length remains finite, consistent 
with a suppression of the phase transition due to disorder. In this paper, we revisit 
this question by combining, in a single experiment, high resolution isotherms and 
light scattering measurements for ^He confined in aerogels close to Tc- 


2 Experimental 

We have studied two aerogels samples with the same porosity (95%) and differ¬ 
ent microscopic structures, obtained respectively by synthesis in basic (B102) or 
neutral (N102) conditions. Both samples are about 3.7 mm thick and 14 mm in 
diameter and located in an optical cell, allowing light scattering measurements at 
45° and 135° from the incident direction-. The cell temperature is measured by 
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Fig. 1 (Color online) Isotherms in B102 close to T,.. The solid line (red online) is the adsorption 
branch, while the dashed line (green online) is the desorption branch. The dashed line in the 
lower right corner (blue online) represents the bulk isotherm— for the given temperature. Bulk 
gas densities are out of range on the scale of this figure and fall between 52 g/L (5.15 K) and 
57 g/L (5.18 K). Hysteresis closes between 5.16 K and 5.17 K. 


a germanium resistor, and regulated to about 10 by a heater glued onto the 
cell walls. Condensation is controlled by changing between 80 K and 200 K the 
temperature of an external reservoir connected to the cell—. The mass condensed 
into the aerogel is precisely obtained from the total amount of "^He in the system 
and the reservoir temperature-. The pressure is measured by a room temperature 
sensor connected to the cell by a separate capillary. During condensation, latent 
heat is released and flows to the cell walls, causing a local increase of the "^He 
temperature, due to its poor thermal conductivity. In order to suppress this thermal 
gradient, we stop condensation at regular time intervals and let the system ther¬ 
mally relax about one hour while continuously measuring the pressure P and the 
heating power W), necessary to regulate the cell temperature. 


3 Isotherms 

Figure[T]shows isotherms for B102 close to 7’c=5.195 K, for condensation or evap¬ 
oration rates between 0.25 and 0.50 ;Umol/s. Each stop is followed by a relaxation 
of the pressure. Once the pressure has relaxed, the loop remains open at 5.15 K 
and 5.16 K while it is closed at 5.17 K and 5.18 K within our pressure resolution 
(10 jubar). Hence, the hysteresis loop closes between 25 and 35 mK below Tc, 
in agreement with the results of Herman et al. on a similar aerogel—. This is in- 
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termediate between the temperatures (respectively 9 and 65 mK below Tc) where 
the bulk correlation length , and the interface width (3.6 ^— ) become equal 
to the estimated aerogel correlation length of 10 nm~ . For the two temperatures 
between the closure temperature Tcu and Tc, the isotherm does not present any 
pressure plateau. As Flerman et al. we do dot find any evidence for any equilib¬ 
rium long range order below Tc. 

The isotherms present a kink at a certain density, above which the pressure 
rises faster, with a slope comparable to that of the bulk isotherms just above the 
saturated pressure, showing that the kink marks the end of the condensation inside 
the aerogel. The corresponding density of the confined liquid is only slightly larger 
than the bulk density (at most 10%). While this agrees with ref.^ and the fact that 
most of ^He is far enough from the silica strands not to be compressed by the 
Van der Waals interaction, this contrasts with the conclusion of ref.—, based on 
neutrons transmission measurements, that the average density of CO 2 in aerogel 
close to the critical point is 50% larger than in the bulk. We will come back to this 
point in 0 


4 Pressure relaxation and slow dynamics 

We have searched for slow dynamics characteristic of disordered systems by ex¬ 
amining in detail the relaxation of pressure following each condensation step. At 
4.8 K, far below Tc, figure |2a) shows that once the condensation is stopped, P 
and W/, relax to their final values with a similar time scale of about 5 minutes. 
This is comparable to the computed decay time (2 minutes) of the slowest mode 
of thermal diffusion in a 3.7 mm thick slab filled with bulk liquid, suggesting that 
the observed pressure relaxation could be related to the cooling of ^He following 
the stop of condensation. Indeed, we have checked that decreasing the cell tem¬ 
perature at constant filling on the pressure quasi-plateau does shift the measured 
pressure by the corresponding change in the saturated vapor pressure (about 1.3 
mbar per mK). 

Figures |2jb) and (c) show relaxations close to Tc, below and above Tch- Wh 
relaxes on the same time scale of 5 minutes with an amplitude proportional to the 
initial flow rate, consistent with a thermal effect similar to that at 4.8 K. In contrast, 
the pressure relaxes on a longer time scale, which depends on the average density 
and is maximal close to full condensation. We can fit the observed relaxation by 
a sum of two exponentials, one proportional to W), due to thermal relaxation, the 
other with a slower decay rate (about 8 minutes at 5.15 K and 12 minutes at 5.18 K 
close to full condensation) which could reflect slow dynamics of the fluid distri¬ 
bution inside the aerogel. This would be consistent with the existence of energy 
barriers. The longer relaxation time at 5.18 K could result from the fact that un- 
surmountable barriers below Tch can be overcome on the observation time scale. 
However, on the one hour timescale of our measurements, we do not observe the 
stretched exponential behavior expected for a glassy dynamics associated with a 
distribution of barriers. Direct measurements of the fluid dynamics using photon 
correlation spectroscopy would be needed to elucidate whether this feature could 
be due to our limited pressure resolution. 
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Fig. 2 (Color online) Relaxations of the pressure P(+) and of the heater power Vf), (x) following 
a condensation step at three temperatures. The inset shows P and Wi, at 4.8 K with the flow rate 
on or ojf. The behavior at 4.8 K is fitted with single exponentials and is attributed to thermal 
relaxation. For 5.15 K and 5.18 K, the dashed and dotted lines are exponentials with the decay 
rates measured at 4.8 K for P and Wh- Pressure relaxation of Wh at 5.15 K and 5.18 K is slower 
than at 4.8 K and can be fitted (solid lines) by the pressure relaxation at 4.8 K, rescaled by the 
amplitude of Wi,, plus a slower exponential, which we attribute to slow dynamics of the fluid 
distribution near 7).;,. 


5 Light scattering close to Tc 

Figure|2a) shows the intensity scattered at 45° as a function of pressure for B102. 
For T < Tch, this signal (as well as that at 135°) increases in the region of the hys¬ 
teresis loop. In N102, the light scattered by the entrance region of the laser sheet 
at 135° shows a similar behavior (the light scattered by other regions, or at 45°, is 
dominated by the strong silica background of N 102^). As at lower temperatures^, 
this behavior is consistent with the progressive growth of microscopic static liquid 
domains as the pressure increases. The decrease of the signal with temperature is 
consistent with the decrease in the optical contrast between liquid and gas as the 
bulk critical temperature is approached. In this interpretation, one could expect 
the optical signal to vanish at 7),. However, measurements in N102 performed 3 
and 50 mK above Tc show a continuous increase of the scattered intensity up to 
the region of maximal slope in the bulk pressure-density isotherm, followed by an 
abrupt drop to the background silica level (Figure [2^b)). Correlatively, the mea¬ 
sured transmission drops in the region of the peak in the scattered intensity (the 
logarithm of the extra attenuation with respect to the bare silica being proportional 
to the scattered signal at 135°, as expected). While this behavior could evocate 
critical opalescence, it can be accounted for by the concept of critical adsorp- 
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Fig. 3 (Color online) Pressure dependence of the light intensity scattered at 45° for B102 (a) 
and 135° for N102 (b), normalized by the intensity Isi scattered by the bare aerogel. In figure (b), 
the curve at 5.25 K is vertically shifted by 0.3 for clarity. Figure (b) also shows the theoretical 
dependence of the excess mass adsorbed on a flat silica substrate, rescaled for comparison to the 
optical signal (see text). 


tion— . The van der Waals attraction by the silica, combined with the large bulk 
compressibility of ^He close to Tc, increases the density of ^He close to the silica 
strands, hence the scattered signal, in a similar way to the case of an adsorbed 
film~. With respect to an homogeneous situation, the scattered field is increased 
by a quantity proportional to the excess mass with respect to the bulk, defined as 
the volume averaged difference between the local density and the bulk density at 
the external pressure. This quantity can be computed as a function of temperature 
in the case of a planar substrate, using the bulk state equation of and the 

van der Waals attraction given in ref.—. As shown by figure [^b), the temperature 
dependence of the obtained excess mass reproduces quite well the observed scat¬ 
tered intensity, the scaling factor being close to that computed from the known 
dielectric constants of silica and ^He, and a typical diameter strand of 2 nm. On 
the other hand, the measured extra scattering due to "^He (about 25% of the sil¬ 
ica contribution, corresponding to a mean free path of 10 mm) is nearly thousand 
times larger than the (separately measured) bulk scattering signal at 7). and 45°. 
Here also, photon correlation spectroscopy would be needed to separate possible 
critical thermodynamic temporal fluctuations from the strong static contribution 
of critical adsorption. 

Finally, we note that our results for N102 are quite similar to those obtained on 
CO 2 using neutrons rather than light scattering—. This suggests that, in the latter 
case, the decrease in transmission measured near Tc close to the critical pressure is 
also due to the increase of scattering due to critical adsorption. Extracting the fluid 
density from the neutrons transmission, as the authors of ref.— do, would then be 
incorrect. This might be the origin of the reported large increase of the confined 
density with respect to the bulk fluid for CO 2 (see 


6 Conclusions 

In summary, in contrast to refs. ^'^^'^^ , our experimental results do not show any 
equilibrium first order transition above 7),/,, the closure temperature of the hys- 
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teresis loop. Either, the transition predicted by the RFIM lies inside the hysteretic 
region below T^h, or the temperature above which the bulk ^He correlation length 
is large enough for the RFIM to apply is larger than T*, the critical temperature of 
the RFIM model. On the other hand, around T^h, the relaxation of pressure does 
not show any clear clear evidence for marked glassy dynamics. Finally, close and 
above the bulk critical temperature, the optical signal is dominated by the static 
contribution resulting from the large compressibility of ^He, masking any possi¬ 
ble critical dynamic contribution. To go further, we plan to improve the sensitivity 
of our measurements by using photon correlation spectroscopy to directly detect 
small dynamic changes of the fluid distribution. 
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